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Abstract. By employing a local two-fluid theory, we investigate an obliquely 
propagating electromagnetic instability in the lower hybrid frequency range 
driven by cross-field current or relative drifts between electrons and ions. The 
theory self-consistently takes into account local cross-field current and ac- 
companying pressure gradients. It is found that the instability is caused by 
reactive coupling between the backward propagating whistler (fast) waves 
in the moving electron frame and the forward propagating sound (slow) waves 
in the ion frame when the relative drifts are large. The unstable waves we 
consider propagate obliquely to the unperturbed magnetic field and have mixed 
polarization with significant electromagnetic components. A physical picture 
of the instability emerges in the limit of large wavenumber characteristic of 
the local approximation. The primary positive feedback mechanism is based 
on reinforcement of initial electron density perturbations by compression of 
electron fluid via induced Lorentz force. The resultant waves are qualitatively 
consistent with the measured electromagnetic fluctuations in reconnecting 
current sheet in a laboratory plasma. 
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1. Introduction 

Current-driven instabilities with frequencies higher than ion cyclotron frequency {cu > 
fli) or wavelengths shorter than ion skin depth {kXi > 1; Aj = c/ujpi) have been a popular 
subject for space and laboratory plasma research [see e.g. Gary, 1993]. Recently, this 
topic has been revisited in the context of magnetic reconnection [see e.g. Biskamp, 2000], 
where intense current density exists locally in the diffusion region. In particular, the Lower 
Hybrid Drift Instability [Krall and Liewer, 1971] (LHDI) driven by a density gradient has 
received considerable attention as a potential source of anomalous resistivity. 

When the LHDI propagates nearly perpendicular to the magnetic field it is purely 
electrotatic. Such waves have been observed at the low-/3 edge of the current sheet in the 
laboratory [Carter et al., 2002a], in numerical simulations [see e.g. Scholer et al., 2003], 
and in space [Shinohara et al., 1998; Bale et al., 2002]. They are driven unstable by inverse 
Landau damping of the drifting electrons. 

However, these electrostatic modes are largely stabilized [Davidson et al., 1977] inside 
the high-/? reconnection layer, where the magnetic field gradient is large and the VS drift 
of the electrons is in the wrong direction to amplify the waves. Further, it is observed that 
their amplitudes do not correlate with the fast reconnection in the Magnetic Reconnection 
Experiment or MRX [Carter et al, 2002b]. By contrast, magnetic fluctuations up to the 
lower hybrid frequency range have been more recently detected [Ji et al., 2004] in this high- 
P center of the current sheet in the MRX. These propagate obhquely to the magnetic field, 
and their amplitudes exhibit positive correlations with fast reconnection. A theoretical 
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explanation for the origin of these magnetic fluctuations, other than the electrostatic 
perpendicularly propagating LHDI waves, is therefore in order. 

Earlier, motivated by observations of high frequency magnetic fluctuations in a magnetic 
shock experiment, Ross attempted [Ross, 1970] the first theoretical exploration of such 
candidate obliquely propagating electromagnetic high-frequency waves driven by a relative 
drift between electrons and ions associated with local currents. Based on a two-fluid 
formalism in the electron frame, Ross showed that unstable waves propagating obliquely 
to the magnetic field are excited by reactive coupling between ion beam and whistler waves. 
Such an instability is generally known as the Modified Two Stream Instability [McBride 
et ai, 1972; Seiler et ai, 1976] (MTSI) since it is driven by a local current across a 
magnetic field unrelated to a diagmagnetic drift. 

Extensions to a full kinetic treatment of both ions and electrons were made for this 
instability [Lemons and Gary, 1977; Wu et ai, 1983; Tsai et ai, 1984]. Unlike the per- 
pendicular LHDI, the obliquely propagating MTSI persists in high-/3 plasmas, where the 
critical values of relative drift for the instability are typically a few times the local Alfven 
velocity, and possesses significant electromagnetic components. However, in most of these 
works, a finite pressure gradient self-consistent with the cross-field current was left out in 
the wave dynamics. This neglect throws doubt on the applicabihty of the MTSI to the 
MRX, where all the current is due to inhomogeneities. 

Recently, global eigenmode analyses [Daughton, 1999; Yoon et al, 2002; Daughton, 
2003] of the current driven instabilities have been carried out to take into account the 
effects of boundary conditions of a Harris current sheet [Harris, 1962]. This followed ear- 
lier work on the same subject [Huba et ai, 1980]. It was found that for short wavelengths 
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{kXe ~ l;Ae = c/ujpe), the unstable modes concentrate at the low-/? edge, and they are 
predominantly electrostatic similar to the perpendicular propagating LHDI. In contrast, 
for relatively longer wavelengths {ky/XgXi ~ 1), unstable modes with significant electro- 
magnetic components develop in the center region. These are similar to the MTSI at 
high-/?. For even longer wavelengths {kXi ~ 1), a drift kink instability [Daughton, 1999] 
is known to exist but this has a slower growth rate at more realistic ion-electron mass 
ratios. More recently, these analyses have been further extended to non-Harris current 
sheets [Yoon and Lui, 2004; Sitnov et al., 2004]. When relative drift between electrons 
and ions is enhanced, the central region is clearly dominated by instabilities resembling 
the MTSI. 

The first numerical simulations of the MTSI have been carried out in a two-dimensional 
local model [Winske et al., 1985], but focused on the electron heating. Particle simulations 
have also been carried out in three-dimensions to study stability of a Harris current sheet 
under various but limited conditions [Horiuchi and Sato, 1999; Lapenta and Brackbill, 
2002; Daughton, 2003; Scholer et al., 2003; Shinohara and Fujimoto, 2004; Ricci et al, 
2004]. It was found that at first the LHDI like instabilities are active only at the low-/? 
edge, and modify the current profile which then leads to the long wavelength electromag- 
netic modes, such as drift kink instabilities or Kelvin- Helmholtz instabilities [Lapenta and 
Brackbill, 2002]. Recent simulations using more realistic parameters (larger mass ratios 
with more particles) in larger dimensions indicate [Ricci et al, 2004] that the MTSI-hke 
modes also develop in the central region. While the characteristics of the observed waves 
in the MRX current sheet are generally consistent with these linear stabihty analyses 
and nonlinear simulation results, there has been yet no convincing physical explanation 
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of the observed electromagnetic waves in the lower hybrid frequency range. Comparisons 
between MTSI and LHDI, the latter of which involves a self-consistent pressure gradient, 
were made based on local kinetic theories [Hsia et ai, 1979; Yoon et ai, 1994; Silveira 
et ai, 2002], but with a focus on nearly perpendicularly propagating waves. Extensions 
to larger propagation angles were also attempted earlier [Zhou et ai, 1983; Zhou and Cao, 
1991] but with few discussions on the underlying physics. 

Motivated by the observations in the MRX and these recent theoretical developments, 
we investigate this instability based on a local two-fluid formalism in this paper. Our 
analysis is of the MRX and includes the self-consistent pressure gradient with large prop- 
agation angles. A local treatment is justified if the wavelength is short {kXi ^ 1) and 
the growth rate is large (7 » Qj), compared to the global eigenmode analyses extending 
throughout the current layer (see for example Kulsrud [1967]). Our focus here is to reveal 
the underlying physics of the instability by using the simplest possible model rather than 
to carry out more involved calculations. We find that when the relative drifts are large, 
the instability is caused by a reactive coupling between the backward propagating whistler 
(fast) waves in the moving electron frame and the forward propagating sound (slow) waves 
in the ion frame. The unstable waves have a mixed electromagnetic character with both 
electrostatic and magnetic components. They propagate obliquely to the unperturbed 
magnetic field. The primary positive feedback mechanism for the instability is identified 
as reinforcement of initial electron density perturbations by an induced Lorentz force. 
The role this instability plays in magnetic reconnection, such as anomalous resistivity and 
heating, will be discussed in a forthcoming paper [Kulsrud et al., 2004] that is based on 
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quasi- linear theory (see also Winske et al. [1985]; Basu and Coppi [1992]; Yoon and Lui 
[1993].) 

2. Theoretical Model 

The basic features of our model are described in this section. Since our main objective 
is to understand physics of the underlying instability, we develop a theoretical model, 
which contains the essential ingredients for the instability, yet remains simple enough so 
that the feedback mechanism can be understood. In contrast to the past work, most of 
which is based on full kinetic theory, we find that we are able to use a simple two-fluid 
theory and still obtain reliable results. We show that most features of the instability can 
be revealed by this simple model. 

2.1. Method of the calculation 

We wish to treat the LHDI mode by an approach somewhat different from earlier 
approaches. Our basic assumption is that the drift velocity is produced by equilibrium 
gradients (LHDI) rather than an ion beam (MTSI). In the MRX the gradients are the 
origin of the relative drift velocity of the ions and electrons which is just the diamagnetic 
currents, so that the instabihty is an LHDI. However, since the LHDI has been usually 
treated as a nearly perpendicular propagating mode and we restrict ourselves in this paper 
to propagation at angles finitely different from 90 degrees, we refer to our instability as 
the obliquely propagating LHDI or more briefly the obhque LHDI. 

The reason we do not consider the LHDI near 90 degrees is that it has been shown 
to be stable in the central regions of the MRX and, as discussed in the introduction, we 
are interested in explaining the observed instabilities there. In fact, as shown by Carter 
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et al. [2002b], the gradient of the magnetic field is large there and the V-B drifts cause 
the resonant electrons to drift in the opposite direction than inferred from their current. 
The oblique instability we investigate is a non resonant one. 

We assume that the mode is at a large frequency compared to the ion cyclotron fre- 
quency, and the wave length is small compared to the ion gyration radius, so that the ions 
may be considered to be unmagnetized. We also assume that the frequency is small com- 
pared to the electron cyclotron frequency, Oe, and that the wave length is large compared 
to the electron gyration radius, so that the electron can be treated by the drift kinetic 
theory. This theory is described in Appendix Al, but the upshot of it is that one expands 
the Vlasov equation in the small parameter pe/ A where A is the perpendicular scale of the 
perturbation as well as the equilibrium. One solves the Vlasov equation to lowest order 
to obtain the zero order electron distribution function /o from which one can obtain the 
electron pressure tensor, Pg. Then one calculates the perpendicular velocity moment of 
the first order distribution function /i, to find the perpendicular electron current. But 
this calculation is equivalent to taking the perpendicular electron fluid equation of mo- 
tion with this pressure tensor. The parallel current is then obtained from the continuity 
equation, V • — d{neie)/dt — 0. 

This procedure is totally equivalent to previous calculations giving identical results in 
the small pe Umit. It might be argued that one should consider waves with k_iPe ~ 1 
since in previous work on the perpendicular LHDI the maximum growth occurs when 
k±Pe ~ 1- However, for the oblique LHDI the maximum growth actually occurs when 
k±Pe <S 1 and the mode becomes stable for k±Pe that approaches unity. (The guiding 
center treatment is appropriate for inhomogeneous systems, since it makes no assumption 
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about near homogeneity and avoids the comphcating approximations concerning it that 
are usually made.) 

It turns out that it is not appropriate to treat the pressure tensor as anisotropic for the 
MRX experiments, in which the magnetic fluctuations are observed. This is because the 
electron-ion collision rate is comparable to the frequencies and growth rates of the mode, 
so it is just as accurate to take the pressure as isotropic. Further, it is also appropriate 
to assume that the plasma is isothermal, so that p — nT in general; po = no{x)T in 
the equilibrium and pi — riiT is the perturbation. The fact that T is constant in the 
equilibrium over the region occupied by the mode is supported directly from observations. 
The fact that perturbations in the temperature are zero follows from the very large thermal 
conductivity along the lines, so that the thermal relaxation time is shorter than the 
perturbation growth time. With this assumption we can avoid the solution for /o and 
work entirely from the electron fluid equation, to determine the perpendicular electron 
currents. 

At this point we are in a position to solve for the ion and electron currents in terms 
of the electric fields. However, one further physical result, charge neutrality, allows us 
to further shorten the calculation. Since the Debye length is very small compared to 
even the electron gyration radius, we may assume to an excellent approximation that 
the perturbed electron density riei is equal to the perturbed ion density rin and this 
enable us to easily evaluate the relevant terms in the perturbed equation of motions of 
the electrons. (Of course if we had avoided this step and solved directly for the ion 
and electron currents separately and then substituted in Maxwell's equations, charge 
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neutrality would have followed automatically. Introducing charge neutrality earlier leads 
to considerable simplicity in the calculation and more physical insight.) 

To summarize our calculation: we first write down the equilibrium conditions. Then 
next we calculate the perturbed ion current and density from the unmagnetized ion dy- 
namics. We then calculate the perturbed perpendicular electron current from the per- 
pendicular equation of motion for the electrons. We can then find the parallel electron 
current from V-j = V - = 0. Knowing these currents, we then substitute them 

into Maxwell's equations to find three independent relations for the wave electric fields. 
However, it turns out that one of the three Maxwell's equations can be simplified to the 
electron force balance along the field line. Thus, this eliminates the needs to calculate 
the parallel electron current directly from V-j = 0, which is demanded by the charge 
neutrality condition. 

2.2. Equilibrium 

For definiteness, we assume the MRX equilibrium is a Harris equilibrium and study it 
in the ion frame. This seems the most physical frame in which to study the instability 
since it turns out to be essentially an unstable sound mode which is carreid by the ions. 
We concentrate our attention on a small region say about half way out from the center of 

the Harris sheet. 

In this frame as shown in Fig. 1(a), there is an electric field Eq balancing their pressure 
force, Tiduo/dy, in the y direction: 

enoEo = T^^. (1) 
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The magnetic field, Bq, is chosen in the z direction. A current is carried by electrons 
drifting in the x direction with a speed Vq. Force balance of the electron fluid then is 
given by 

-eno{Eo - VoBo) = (2) 
Eliminating drio/dy in Eqs.(l,2), we have 

Eo = -^VoBo. (3) 

-'■e ' -'■i 

If the plasma resistivity is finite, the electron current in the x direction cannot be main- 
tained without an electric field in the same direction, E^o- We shall see later, however, 
that its effects on the wave dynamics are small as in the MRX. 

3. Dispersion Relation 

All wave quantities are assumed to have a normal mode decomposition proportional to 

exp[i(k • X — cot)] 

with the wave vector k = {k^, 0, k^) and the wave angular frequency uj. Note that k here 
does not have a y component. This assumption is justified in a local theory if wavelengths 
are much smaller than the current layer thickness in the y direction. 

The governing equation between cu and k, or the dispersion relation, follows from three 
independent equations that relate the three components of the wave electric field, E^, Ey, 
and Ez- These can be derived from Ampere's law and Faraday's law, 

k X (k X E) = -icufxoi, (4) 

which leads to 

k^E^ - k^k^E^ = iuJUojx (5) 
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k'^Ey = iuJiiojy (6) 
klE^ - k^k^E^ = iuJI^ojz- (7) 

Here /iq is the vacuum magnetic permeability. Next, we separately consider ion and 
electron dynamics to express the above equations in terms of the electric field. 

3.1. Ion Dynamics 

We take the ions as unmagnetized and solve the kinetic equation for the perturbed dis- 
tribution function assuming the equilibrium ion distribution function is Maxwellian with 
constant temperature, but variable density, drio/dy = cuq. From Eq.(l), e = eEg/Ti = 
2eEo/Mvf. 

The solution of the ion Vlasov equation is carried out as an expansion to first order in 
e. The result is most easily expressed in terms of the electric field components Ei and E^ 
defined in Fig. 1(b), in which Ei is the component parallel to k, and £'3 is the component 
perpendicular to it and in the x — z plane. The perturbed ion current can then be written 
(Appendix Bl), 

- i(e/k)((Z' + Z)E,k\ (8) 

and the perturbed ion density is 

r^ = ^^fiOi^-^ + ieEy) (9) 

where k = k/A;, C, — u/kvi, and Z is the plasma dispersion function. We find that for the 
principal instabilities the phase velocity is somewhat larger than Vi so for convenience we 
first take the C ^ 1 limit (the cold limit), determine the parameter range of instability. 
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Then, in Appendix B3, we are able to employ a simple modification of the dispersion 
relation to extract the correct growth rate including the finite ion thermal effects. 

In the cold limit the ion current neglecting the e correction is obtained from the C ^ 1 
limit and is 

r ^ i^eoE, (10) 



where the ion plasma angular frequency Upi = yrioe^/Meo and eo is the vacuum suscep- 
tibility. In the same limit the perturbed ion density is 



The neglected e term is much smaller than the other one since, for our local theory, we 
assume k/e ^ 1. Indeed, it is shown in Appendix B2 that the neglected term only has a 
small effect on the dispersion relation. 

3.2. Electron Dynamics 

As we have shown in Appendix Al, the perpendicular electron current can be obtained 
from the first order force balance for the electron fluid, 

j"" X Bo = enoVo x B + cuqE + enEo + TgVn + mriQ— — (12) 

at 

where U^; = E x Bo/i^o ^^'^ ^ electron mass. As shown in Appendix A2, the 

electron inertial terms contribute a small effect to the distpersion relation and we can 
neglect them when determining the instability. The y and x components of Eq.(12), 
therefore, are given by 

-jIBq = -enoVoB^ + enoEy + enEo (13) 
j^Bo = enoE^ + ik^T^n (14) 

DRAFT February 5, 2008, 2:24pm DRAFT 



X - 14 JI, KULSRUD, FOX, AND YAMADA: OBLIQUE ELECTROMAGNETIC DRIFT INSTABILITY 

respectively. Here, = kxEy/uj. Since Eq and n are given already by Eq.(3) and Eq.(ll), 
and jy can be expressed in terms of the electric field. 

We note on the righthand side of Eq.(14) that there would be another term, enE^o, 
where E^o is the unperturbed electric field. However, we will treat it as second order and 
balanced by quasi- linear terms [Kulsrud et ai, 2004]. In fact, the contribution from this 
term is small when compared with the last term if kx ^ eExo/Te as is often satisfied in 



The z-component of the electron current, is not determined by Eq.(12). It turns out, 
however, that it is unnecessary to explicitly calculate it in order to obtain the dispersion 
relation due to simplifications of the 2;-component of Maxwell's equation, Eq.(7). This is 
because the electrons are so easily accelerated along the field line by the force, F^, on the 
electron fluid where 



The various terms in this force are separately large and must balance closely to avoid very 
large parallel electrons currents. In fact, taking j| = —noev^ = —i{e/mu!)F^ and using 
Eq.(lO) for the ion current, we can write the z-component of Maxwell's equation, Eq.(7), 

as 



Since uj^^/uj^^ — m/M <C 1 and {k\e)^ -C 1 to our interests here, the above equation 
simplifies to the one demanding the electron force balance in the z direction. 



the MRX. 




klE^ - kxk^Ex = iujfJ'oijl + jl) 




(15) 



E. + VaBy + ik, — — ^Q, 



e no 



(16) 
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where By can be expressed in terms of the electric field using Faraday's law, 



B,. — 



UJ 



In Appendix A2, we show that the neglected terms have only a small efi^ect on the dis- 
persion relation. We note that, although unneeded for the dispersion relation, the z- 
component of the electron current, can be determined by V • j = 0. This is a conse- 
quence of the charge- neutrality condition, which is in turn enforced by Eq.(16). 

It is interesting to note that if we allow the propagation angle approach to 90°, the 
parallel phase velocity can be comparable to the electron thermal velocity. In this case, 
we need to include a Landau term in Eq.(16). Then, if /3e <^ 1, we would be able to 
recover the electrostatic perpendicular LHDI [Krall and Liewer, 1971]. However, since 
this electrostatic LHDI disappears at the high-/? of interest to us, we need not include the 
Landau term. 

3.3. Dispersion Relation 

Substituting expressions of j^, j\ and n [Eqs.(10,13,14,ll)] into Eqs. (5,6,16) we obtain, 
after some algebra, the dispersion relation 



^ -^xx -^xy -^xz ^ 

-Dyx -^yy -^yz 
\ ^zx ^yz ^zz J 

where 

Dxy^ i(n-KV sine) 

r.9 ■ n n A KV COS 6 

Dx,= -K^ sine cose --^ — - — 
Dyx- -M^-y^^ 



f Ex\ 

Ey =0, (17) 
\Ez) 
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D 



yy 



D 



.(3eK^ sine cose 



yz 



KVcose 



2 Q. 

(3eK^ sin ^ COS ^ 



D 



zy 



n-KVsme- 



/?e K'^ COS^ e 



Here the dimensionless parameters are defined by 



a^^, K^k^, i/.ii, A. "o^- 



UJ, 



pi 



Pi = p2 /o ' sin^ = — . 



[18) 



Here, cUd is the ion cyclotron angular frequency cBq/M and Va the Alfven speed 
So/^//^Mr^^. 

The KV term in D^x and in D^z and the (3e terms all result from replacing the kinetic 
equation for the perturbed density n by it cold limit. The 'one's in 

Dxx and Dyy are ion 

currents which are similarly appoximated. 

The resultant dispersion relation Q,{K) is a fourth order algebraic equation in Q with 4 
controlling parameters, V, (3e, Pi, and e, 



{K^ + 1) {K^ cos' e + 1)- K'V' sin' e + ^K' 



+KVsine 
+K' 



n 



2 



/9e + A 

{K' + lfcos'e-K'V'sin'e 



-{K' + l)V' 



/5e + A 



0. 



(19) 
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4. Wave Characteristics and Instability 

4.1. Basic Wave Characteristics without Drift 

The basic wave characteristics described by Eq.(19) are summarized here for the case 
that there is no drift between ions and electrons. When V — and ^ = 0, Eq.(19) reduces 

to 

which represents four waves, as shown in Fig. 2 for the case oi Pe — Pi — ^- Two waves 
are whistler waves, traditionally termed fast waves, while the other two waves are sound 
waves or slow waves. One of each waves propagates along the background magnetic field 
and the other propagates against. As expected, the whistler waves are largely transverse 
waves or electromagnetic waves since the electric field vectors are perpendicular to the 
propagation (k) direction 4> — 90°, where cos0 = k • E/(|k||E|). In contrast, the sound 
waves are largely longitudinal waves or electrostatic waves since ~ 0. 

The situation changes when 9 and /? are varied. In Fig. 3, the angles between E and 
k, 0, are shown ior V — and a few cases of 9 and p. It can be seen that when 9 is 
larger, the whistler waves become less electromagnetic and more electrostatic while the 
sound waves become more electromagnetic and less electrostatic. This trend is stronger 
for larger values of p. 

4.2. An Obhque Electromagnetic InstabiHty 

It is evident that the whistler waves are supported by fast electron dynamics while the 
sound waves are supported by slow ion dynamics. When there is no drift between these 
two fluids, all wave branches stay separate in the dispersion diagram as shown in Fig.2 for 
^ = 0. The situation is similar for more general cases of ^ 7^ 0. liV — 0, Eq.(19) reduces 



Q2 _ 
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to 



{K' + l){K^cos^e + l) + (^K' 



+^K\K^ + If cos" e^Q, (20) 
2 

which represents four waves in the left panels in Fig.4 for the case of 6* = 60° and /5e — 
Pi — 1. It is seen that at this propagation angle, is ~ 40° for whistler waves and ~ 0° 
for sound waves. 

When there is a finite electron drift in the ion rest frame, the whistler waves are doppler- 
shifted so that each fl from Eq.(20) is increased by i^'l^sin6', shown as dotted curves in 
the top-right panel of Fig.4 for the case oiV — 6. In constrast, sound waves, unaffected 
by the drift, are shown as dotted straight lines. When the drift is large, some part of the 
backward propagating whistler waves branch can intercept with the forward propagating 
sound wave branch, resulting in instabilities through reactive couplings. The case oiV — 6 
is shown in the right panels of Fig.4 and all other parameters are the same as in the left 
panels. It is seen that when K <~ 6 or X >~ 16, all four roots are real and thus all 
waves are stable. When 6 < < 16, two of roots become complex conjugates as a result 
of coupling; one of them is damped and another growing (the growth rates are shown in 
the middle-right panel). The maximum growth rate is about 8 times of a;cj at ~ 11. 
Since the polarization angle (f) ~ 15°, the unstable waves have significant electromagnetic 
components. 

Figure 5 shows the unstable region and contours of polarization angle in the 9 — K plane 
for a few values of V. It is seen that the unstable waves are localized to small K when 9 
is small and to large K when 9 is large. The unstable region expands and the growth rate 
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increases with increasing V. The polarization angle ranges between 10° to 25°, and is 
larger near the small K and small 9 corner. 
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5. A Physical Picture 

5.1. Further Simplification of Electron Dynamics 

In order to understand the primary feedback mechanism of our instability we make 
further simphfications to the dispersion relation given by Eq.(17). We first start by 
rotating the coordinate for E as shown in Fig. 1(b): {E^, Ey, Ez) to {Ei, E2, E^). Ei 
is in the k direction, representing the electrostatic component. E2 is the same as Ey 
and E3 is another perpendicular component to k, and both of these are electromagnetic 
components. Using the new bases, (Ei, E2, E^), Eq.(17) reduces to 

= 0, (21) 

\ ^31 ^32 ^33 / \Es J 

where 

• n A KV 

Di2^ i{n-KV sine) 
Dis^ -(K'^ + l)cose 
n ^^^^^ (a^ f^^T<'A 

D22 = K'' + l 

D23 — ifi cos 9 

Ds2^ 

D33 = QsmO- KV. 
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Again the KV term in Dn and the /5e terms result form approximating the perturbed 
density, and the 'one's in D12 and D22 from approximating the ion currents. 

Next we simphfy these equations by taking the hmit of large fl, K, and V since this 
asymptotic limit will make the physical mechanism of the instability clear. The simplified 
matrix then reduces to 

/ -iKV sine -K-'cose^ 



n \ 2 J 



(22) 



\ n \ 2 J J 

Each line of the above matrix equation represents the balance of the leading forces on 
the electron fluid along the three coordinate directions y, x, z, respectively. By referring 

back to Eq.(7) and Eq.(12) we can see that the force balance can be written 

y : -enEo - jo^B^ - J^Bq = 

X : —erioEi sin 9 — dpe/ dx + jyBo = , (23) 
z : —erioEi cos 9 — dp^/dz + joxBy — 

where in the asymptotic limit the current j is all due to the electrons. Interestingly, the 
electrostatic force is balanced by the Lorentz force in all directions. In the y-direction, 
the unperturbed electrostatic field acting on the perturbed electron density is balanced by 
the Lorentz force, which consists of both magnetic pressure gradient, —joxBz, and tension 
—jxBo forces. By contrast, the perturbed electrostatic field is balanced by the magnetic 
tension, joxBy, in the 2;-direction. 

5.2. The Case of 0=0 

We start with the simplest case, ^ = 0, in which there are no perturbed forces in the 
x-direction. In the y-direction the perturbed magnetic pressure force is also zero since 
Bz — kxEy/uj — 0. Therefore, the electrostatic force, —enEg, must be balanced by the 
magnetic tension force, —JxBq. Suppose that the electron density is perturbed in a way 
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such that n > at the origin as illustrated in Fig. 6 (a) in the y — z plane. Because Eq 
points in the positive y-direction, the perturbed electrostatic force on the electron fluid, 
—enEo, points in the negative y-direction at the origin. Since it varies in z, this force 
bends the field line until its magnetic tension force —jxBo balances the —enEo force. 
(Here the field-line bending can also be understood as a result of the perturbed due to 
changing the number of the charged carriers by the perturbed density n.) 

In the 2;-direction, there is now a component of the magnetic tension force towards the 
origin joxBy due to the bent line, as illustrated in Fig.6(a). This force reduces or reverses 
the perturbed electrostatic force —cuqEi produced by the electron density perturbation. 
In the latter case, the perturbed electrostatic force is directed away from the regions where 
n > and towards the regions where n < 0. As a result the perturbed electric field, Ei, 
must point from the regions where n < to the regions where n > 0, such as the origin. 

To see that this leads to instabihty consider the ions which only see the electrostatic field 
El. This electrostatic field will force the ions to condense further at the origin increasing 
their density perturbation. By charge neutrality this will increase the initially assumed 
electron density perturbation and thus lead to instability. 

5.3. The Case of > 

We find that it is convenient to take the limit of Pe — ^ for the discussion of this 
more general and complicated case. Here the feedback to initial perturbations through 
compression or decompression of the electron fluid along the 2;-direction is unaffected 
except for a reduced efficiency. However, there are perturbed forces in the x-direction. As 
before, we suppose an electron density perturbation n > at the origin. When the mode 
is unstable, the perturbed electrostatic force, which is parallel to k, has an x-component. 
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—erioEi sin 9, pointing away from the regions where n > towards the regions where 
n < also as before. This force on the electrons decompresses the magnetic field in the 
n > regions and compresses it in the n < regions. This is illustrated in Fig.6(b) in the 
X — z plane. Because k makes a finite angle to Bq, the magnetic field lines are distorted 
to have both a tension force and also a magnetic pressure force. Therefore, must be 
negative (decompressed) at the origin where n > and thus, the associated magnetic 
pressure force in the y-direction, —joxBz, is directed towards positive y-direction. As a 
result, this force counters the initial electrostatic force, —enEo, (which bends the field 
line) and thus, reduces the tendency towards instability. 

Both these stabilizing and destabilizing forces are included in the dispersion relation 
from Eq.(22), in which we restore Pe to obtain 

Consider a given (large enough) V, it can be seen that instabihty occurs when K exceeds 
some threshold values, and stability returns eventually in the limit of large K, consistent 
with Fig. 4. Thus, if pe is small enough, the growth rate reaches its peak at a wavelength 
longer then pe- However, it is clear from the above equation that, if (3e — 0, the instability 
persists over all K above its critical value (at least until some finite electron gyroradius 
effects become important.) Prom this, we can see that our calculation is essentially based 
on a two-fiuid model, and it is not strictly a Hall MHD calculation since the ions are totally 
unmagnetized and one cannot set /3e = without losing the above physical contents. Our 
calculation is perhaps closer to a hybrid model [see Birn et ai, 2001] with kinetic ions and 
a massless electron fiuid, but in three dimensions. We emphasize here that the background 
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ion pressure gradient is essential for the instability in both 9 = and 9 > cases because 
of the important role played by the associated equilibrium electric field, Eq. 

6. Discussions and Conclusions 

In the MRX, it has been observed that the usual electrostatic LHDI, propagating per- 
pendicularly to the magnetic field, is active only in the low-/? edge of the reconnection 
region, but not in the high-/? central region [Carter et al., 2002a]. This is consistent with 
the theoretical prediction that the perpendicular LHDI is stable at the high-/? [Davidson 
et al., 1977; Carter et al., 2002b]. On the other hand, it has been found that, in the high-/? 
central region, obliquely propagating, electromagnetic waves in a similar frequency range 
are active, and their amplitude positively correlate with the reconnection rate [Ji et al., 
2004]. Motivated by these observations, we have developed a simple two-fiuid formalism 
to derive and analyze in detail an electromagnetic drift instability in the lower-hybrid 
frequency range. We term this the oblique LHDI. 

We show that the main features of the instability are consistent with fully electromag- 
netic kinetic calculations [Lemons and Gary, 1977; Wu et al., 1983; Tsai et al, 1984]. 
We find that, contrary to the perpendicular LHDI, the oblique LHDI persists in high-/? 
plasmas. Further, the growth rate peaks at longer wavelength than electron gyroradius, 
justifying our assumption that the electrons are magnetized. The resultant waves have 
mixed polarization and significant electromagnetic components. The instability is caused 
by reactive coupling between the backward propagating whistler (fast) waves in the mov- 
ing electron frame and the forward propagating sound (slow) waves in the ion frame, and 
occurs when the relative drifts are large. After further simplifications of the model, the 
primary positive feedback mechanism is identified as a reinforcement of initial electron 
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density perturbations by compression of the electron fluid by an induced Lorentz force. 
Interestingly, the revealed mechanism of the instability requires close interactions between 
the electrostatic and electromagnetic forces. In contrast to most of previous theories on 
MTSI, our analysis also suggest that the self-consistent background-ion-pressure gradient 
is essential for the instability. 

A few comments on three-dimensional particle simulations are in order. In addition 
to the dimensionless parameters of Eq.(18), the mass ratio, M/m, is another important 
parameter. To make the simulations feasible, often M/m is hmited to a few hundreds. 
In contrast, our analysis based on the above simple local model is valid in the limit of 
large M/m since ions are treated as unmagnetized. Small mass ratios used in simula- 
tions will limit available wavenumber window for the instability due to the condition of 
A^^ <^ k <^ X~^. In addition, the limited grid size and resolution may not permit nu- 
merical treatment of the large oblique wavenumber range where our instability resides. 
Future numerical simulations with increasingly powerful computers may help to elucidate 
these effects more clearly especially with regard to nonlinear consequences in magnetic 
reconnection. Simulations of non-Harris current sheets, as attempted in the linear analy- 
ses [Yoon and Lui, 2004; Sitnov et al, 2004], may prove to be more physically meaningful 
since they may represent the reality more accurately. 

Many of the predicted features of unstable waves discussed in this paper are also qual- 
itatively consistent with the observed magnetic fluctuations in the MRX [Ji et al, 2004], 
including their existence in the high-/5 region, their frequency range, and their propaga- 
tion direction with respect to the background magnetic field. In fact, the parameters we 
use in the calculation have been drawn directly from the MRX experiments, and they are 
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valid throughout the bulk of the MRX current sheet. Also, the instability does indeed 
persist into the /3e ^ 1 regimes, but the physics of the instability is still uncertain in the 
region where the magnetic field nearly vanishes. One particular comment on their phase 
velocity is worth making. The experimentally measured phase velocity is of the same 
order as the relative drift velocity. Even given the large experimental uncertainties such 
as the measurement location and the unknown relative velocity between the ion frame 
and the laboratory frame, the measured phase velocities are considerably larger than our 
theoretical predictions. As seen in Fig. 4, the unstable waves should have phase velocities 
on the order of the ion thermal speed. However, the theory presented here is limited to 
the case where ky = 0. The phase velocity may be substantially increased by incorpo- 
rating a nonzero ky. This is a subject for future work. Increasing the phase velocity to 
values much larger than ion thermal speed may also help mitigate another shortcoming 
of our analysis: the reduction of the growth rates by ion thermal effect. The role which 
this instability plays in magnetic reconnection, such as in the production of anomalous 
resistivity and its effect on heating, is discussed in a forthcoming paper [Kulsrud et al, 
2004] that is based on quasi-linear theory. 

Appendix A: Detailed Calculations of Electron Dynamics 
Al. Drift Kinetic Equation for Electrons 

Normally, the drift kinetic equation is developed for both electrons and ions, and is 

combined with Maxwell's equations to achieve some important simplifications. This full 

formulation is described in a number of places, for example in Kulsrud [1983]. However, if 

the ions are unmagnetized, as in this paper, the formulation is reduced to that of solving 

the electron Vlasov equation alone, as an expansion in Pe/X, and l/cUcet where A is the 
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length scale of the phenomena, and t is its time scale. We follow the procedure given in 
the handbook article. It is clear that the electronic charge can be used as a guide to the 
expansion and we use 1/e as the expansion parameter. 
The electron Vlasov equation is 

df e 

^ + V • V/ - - (E + V X B) • Vv/ = 0. (Al) 
at m 

We first carry out the expansion for the full distribution, (equilibrium / and perturbed 
5f) and later carry out the expansion in the instability perturbation. 

The lowest order Vlasov equation is accordingly 

-— (E + V X B) • Vv/o = 0. (A2) 
m 

We introduce the E x B velocity by 

Uk = ^Ij? (A3) 

and carry out the transformation of the velocity at each point r, 

V = UE(r) + v' — \Je + v± cos 0x' + v± sin 0y' + v\\h (A4) 

where x', y' and b are local coordinates at each point r, and v±, (j) and v\\ are cylindrical 
coordinates for v'. Then Eq.(A2) becomes 

eBdfo _ eE^dfo ^ ^ 
m d(j) m dv\\ 

If E\\ is non zero, /o would be constant along a helical orbit in velocity space that extends 
to infinity, which is impossible. Thus, E'y must vanish to lowest order and E\\ must be 
considered first order. 

Dropping the second term we see that /o is independent of 4> (gyrotropic) and, thus, a 
function only of t, r, v±, and v\^. 
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Proceeding to next order in 1/e we get 

eB dfi fdfo \ e ^ dfo , . . 

= ^ + V • V/o - -E. — (A6) 



m dcf) \ dt J m " dv\j 

where the expression in parentheses must be transformed to t, r, v±, v\\, (j) coordinates. 

Equation (A6) can only be solved for /i if its average over (which ehminates dfi/dcf) 
) vanishes. The result is 



^ + (Ue + ^;||b)-V/o 

- y(V-U^-b-VU^.b + ^||V-b)|^ 



(AT) 



where DUE/Dt = dlJ E/dt + (U^ + hv\\) ■ We- (Note that the eE\\ term is zero order 
since E\\ is first order and e is minus first order.) 

In principle, /o can be solved for from this equation. For the case of the instability, 
/o can be written as /q + Sfo where /o is a local Maxwellian, is a perturbation, and 
Bo = BqZ to lowest order. The only equilibrium term that survives is the I'yb • V/q term 
so the only restriction on /q is that it be constant along the magnetic field. 

To get the electron current perpendicular to B we need /i, 

•v±fid(f)v_Ldv_Ldv\i = e J (A8) 

which can be obtained directly from Eq.(A6) by multiplying it by dv±/d(j) — —v± sin 0x+ 

v± cos (py, dividing by B and inegrating over velocity space. In fact, we could just as well 

have multiplied Eq.(A6) by and integrated it to find jg x B. Even simpler, we could 

have multiplied Eq.(Al) by B integrated over velocity space and taken the perpendicular 
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part of the result. This result would be the perpendicular part of 



Here the stress tensor is zero order, and can be found from /o once we have solved Eq.(A7) 
for it. 

If we inspect Eq.(A9) we see that the inertia term and the V • P are zeroth order, 
but the neE term is minus first order in the 1/e expansion. Thus, has a minus first 
order part, ne time the E x B drift, and zero order parts, essentially the diamagnetic and 
polarization currents. If the ions were magnetized, this minus first order current would 
be cancelled by the corresponding E x B current of the ions, but this is no longer the case 
for unmagnetized ions. 

This procedure gives the perpendicular current of the electrons. The parallel current is 
given by the continuity condition 



Again for finite n the ne term is minus first order, no is given by the zero moment of 
/o. However, ui is needed to give the finite parallel electron current, and for it we need 
the zero moment of fi. This zero moment cannot be obtained from Eq.(A6), which only 
gives the dependent part of /i, dfi/dcp. To get the mean part it is necessary to go to 
next order in the 1/e expansion of the Vlasov equation. This has been done some time 
ago [Frieman et ai, 1966], and will yield ni. 

This procedure is certainly possible to carry out in all detail as outlined above and is 
fairly easy for our perturbation problem. In fact if it is carried out in a velocity frame in 
which the equilibrium electric field is zero (the so-called Harris frame) the results turn out 
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to be essentially identical to those calculated by Yoon et al. [1994] in the common limit 
of approximation, small gyration radius and small frequency compared to the electron 
cyclotron frequency. As stated in the text, we can avoid some of the calculation by taking 
the perturbed density from that of the ions by quasi neutrality. This also avoids going 
to next order in the Vlasov equation to find rii. This assumption puts a constraint on 
E\\ in an early phase in the calculation rather than waiting for substitution in Maxwell's 
equations to enforce it. In any event the drift kinetic approach is completely consistent 
with earher calculations of LHDI. 

A2. Electron Inertial Terms 

The first two rows of the matrix in Eq.(17) represent —ft/rioe times the y and x com- 
ponents of Eq.(12). Their initial terms are i{m/M)noeflEx and —i{m/M)noeflEy, re- 
spectively. Multiplying these by (— Q/noe) we get for the first two rows of the matrix 
equation 



DyX Dyy + Dyz 



where the coefficients D are given by Eq.(17) as before. The last row represents —Qc^/cOp^ 
times the last term in Eq.(15). Bringing all the other terms to the right-hand side and 
multiplying these by -c^/u^^, we get -{m/M)K'^ sin 9 cos OE^ + (m/M) (1 + K"^ sin^ e)E^. 
Multiplying these by VL and adding the results to the last row of the matrix equation, we 
obtain 

TYl Ttl 

K^flsinecose D,y D + —^(1 + sin'^ 9) . 
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Transforming to the {Ei, E2, E^) components of the electric field, we have 

TJX Tfl 

Du-^^^sm9 D12 Di3 + i—n^cos9 

M m ^ 

D21 D22 + —fl^ D23 

777 ?77 

Dsi + —ncos9 D32 Ds, + —nsme{l + K^) 
and for the limit of large K, V and Eq.(22) becomes 

^ --T^^-Z^'sine -tKVsine - cos 6 + i^n^ cos 6 ^ 
/3e + f3i Q M M 

_^sin^/,_^ \ + 

n \ 2 I M 

^(n2_^^2V^f^cos^ -KV+'^nsineil + K^) 

\ n \ 2 J M M ^ ' J 

If we regard K, V, and ^2 as all of order K, then we can we can see that the relative 
corrections are of order at most m/M except in the one-one and three-three elements 
where they are of order ~ {m/M)K. These corrections are all small and can be neglected 
as long as X ^ M/m. 

Incidentally, the correction in the third line represents the extra parallel electron field 
needed to accelerate the electrons along the magnetic field to achieve charge neutrahty. Its 
smallness indicates the ease with which the electrons are able to achieve charge neutrality. 

Appendix B: Detailed Calculations of Ion Dynamics 
Bl. Perturbed Ion Current and Density 

The expressions for the unmagnetized ion current and density given in Eqs. (8) and (9), 

which keep the equilibrium density gradient, as a first order correction are found from the 

perturbed ion distribution function with the same correction. The latter is obtained by 

iterating the perturbed ion Vlasov equation 

- k • v)A + + -E,— + _E, . — = (Bl) 
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The second and third terms are the correction terms. Therefore, drop them at first and 
solve for the uncorrected fi from the remaining equation, in the standard way. 



^ -inoe 2vEi e'^^^^ 
■'^ kMvUvz-oj/k) 7r3/2 ^"^^f 

where — 2Tj/M , and where without loss of generality we take the z axis along k. 
We see that dfi/dy = e/i where e — {dno/dy)/no. 

Next, we insert this expression into the second and third terms of the full Vlasov equa- 
tion and solve for the correction, 6f, to /i which satisfies 

- k • v)5/ = -Vyeh - ^Eo^^ (B3) 

After some algebra we can express the zero and first moments of (/i + 5/) in terms of 
the plasma dispersion function, Z, oi — uj/kvi and, thus, obtain Eqs.(8) and (9) of the 
main text. 

B2. Dispersion Relation with the Correction from Background Density 
Gradient 

In Eq.(ll), a term proportional to the density gradient has been neglected in deriving 
the dispersion relation. It is straightforward to show that, by including this term, the 
dispersion matrix is given by 



D D + ^^""^ D 

^ ^ A KVcose ^ 



(B4) 



where the coefficients D are given by Eq.(17). The resultant dispersion relation remains 

as a fourth order equation, and the added new terms only have a small effect on the 

solutions. In the right and middle panel of Fig.4, the growth rate by Eq.(B4) is shown as 
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the dotted hne, which differs iittie from the solid line by Eq.(17) especially in the large K 
limit. (The dotted line indicating instability at very small K has no physical significance 
since the local approximation becomes clearly questionable for such cases.) 

B3. Growth Rates with Warm Ions 

The most important instabilities occur for very local perturbations with large V and 
O. We restrict the discussion of the thermal corrections to this case. 

The cold ion approximation involves using Eq.(ll) for the ion density instead of Eq.(9) 
and Eq.(lO) for the ion currents instead of Eq.(8). In equation (18) the ion currents are 
negligible and only the one-one element and the j3(. terms are proportion to the perturbed 
ion density. Thus the matrix of Eq.(22) with the corrected ion density is 



/ 



-a- 



A 



/3e + A 



-KV 



-isinein^-^K^a 



cos^fo^-^W 



-iKV sine -K^cosO^ 
^2 

-KV 

I 



(B5) 



where 



a = = C ^ (C) 

n 



where C = u;/{kvi) = n/{K^/pl). 
The dispersion relation from Eq.(B5) can thus be written 



_j^4y2^J^^ + kW^ sin^ 9 h^-a^K^ 



Pe + Pi 



K^cos^e(n^-a(^KA =0. 



(B6) 



(B7) 



By dividing this equation by a we see that fi^/a satisfies the same equation as fll, the 
approximate solution for the growth rate with cold ions. Thus we can write 
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where Co = ^0/ (-^^A)- Thus from Eq.(B8) we plot ratio of the true ( to the approximate 
Co as a function of Co in Fig. 7. (Actually, the Q's are pure imaginary so we plot F/ (Ky/Piys 
where the F's refer to approximate and exact normalized growth rates.) 

We see that there is indeed a difference of order unity between the approximate and 
exact values of O or C- Since we see from Fig.4 that the peak Tq/K 1, the true 
r/K ^ 0.5 when Pi — 1 and T/K ki 0.7 when Pi — 0.5. In spite of this reduction, we see 
that the obhque LHDI is still unstable. 
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Figure 1. (a) Illustrations of the equilibrium state. Ions are at rest while electrons 
drift toward positive x direction, crossing magnetic field in the z direction. The resultant 
Lorentz force and electric field is balanced by pressure gradients in the y direction, which 
points towards the current sheet center, (b) Definitions of Ei and E3. E2 is same as Ey. 
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Figure 2. (top) Dispersion relation for the case that V — 9 — and (3e — (3i — l. There 
are 2 whistler (fast) waves (solid lines) and 2 sound (slow) waves (dotted lines), (bottom) 
Angle (0) between E and k vector for both whistler waves (sohd hne) and sound waves 
(dotted line). 
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Figure 3. Angle between E and k for the cases of ^ = 0, 45°, 85° and /3e(= A) = 1 and 
10. Solid lines represent whistler waves and dotted lines represent sound waves. 
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Figure 4. Dispersion relation (top) for the case of no drift (left) and large drift (right). 
Growth rate (middle) and 4> (bottom) are also shown for both cases. See main text for 
the detailed explanation of the dotted hnes. 
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Figure 5. Unstable region where Im(Q) > (filled regions in top panels) and contours 

of polarization angle (0, bottom panels) in the — K plane for the cases of = 3, 6, 10 

and = Pi — 1. 
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Figure 6. Illustrations of the instability mechanism (a) when ^ = in the y — z plane 
and (b) when > in the x — z plane. This distance between the field lines indicates the 
relative field strength. 
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Figure 7. 



Ratio of exact growth rate to approximate growth rate as a function of 



DRAFT 



February 5, 2008, 2:24pm 



DRAFT 



